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Abstract. We investigate the statistical evolution of magnetic neutron 
Q^' stars recycled in Low Mass Binary (LMB) systems, simulating synthetic 

^ I populations. Irrespective to the details of the physical models, we find to 

^ ■ be significant the fraction of neutron stars spinning close to their mass 

I shedding limit relative to the millisecond pulsar population. The sim- 

■ ulated neutron stars show a tail in their period distribution at periods 
shorter than 1.558 ms, the minimum detected so far. Crustal magnetic 
field decay models predict also the existence of massive rapidly spinning 
neutron stars with very low magnetic moment. 



1. Magnetorotational evolution scenario 

Population synthesis models are evolved within a simple recycling scenario (as 
outlined e.g. by Lipunov 1992), where the neutron stars (NSs hereon) may 
experience the phase of ejector, accretor or propeller. The magnetorotational 
evolution includes the general relativistic effects as modelled in Burderi et al. 
(1999). We considered two possibilities bracketing uncertainties in the accretion 
histories of the NSs in LMBs: (i) steady accretion for a time trlo] (h) per- 
sistent accretion followed by a transient phase during which the mass transfer 
becomes unsteady, mimicking the quenching of accretion (Ergma et al. 1998). 
In particular we account for the angular momentum losses by propeller, both 
for persistent and for non-stationary accretion, during the Roche Lobe Overflow 
phase. 

The NSs in our population synthesis models have an initial gravitational 
mass of 1.4 M©. Cook, Shapiro & Teukolsky (1994) showed that all viable equa- 
tions of state (EoSs) for nuclear matter allow for recycling to ultra-short periods. 
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in the case of a nonmagnetic NS. The minimum period below which the mass- 
shedding instabihty hmit is encountered, depends however on the equation of 
state. We consider either a soft (FP) and a stiff (PS) EoS, whose critical periods 
for mass-shedding are Pgh = 0.73 ms (for FP-EoS) and of Pgh = 1-40 msec (for 
PS-EoS), respectively. 

A basic assumption for the present investigation is the crustal origin of 
the NS magnetic field. Its subsequent evolution is governed by the induction 
equation for ohmic diffusion and advection (see Urpin ct al. 1998). While 
at the surface the standard boundary condition for a dipolar field applies, the 
inner boundary condition is subject of scientific debates. Hence we explored two 
hypotheses, imposing at the crust-core boundary either complete field expulsion 
by the superconducting core (Boundary condition I, BC I) or advection and 
freezing (Konar &: Bhattacharya 1997) in a very highly conducting transition 
shell (BC II). 

In paper I, Possenti et al. (1998) evolved NSs in low mass binaries (LMBs) 
through the five phases of the recycling scenario, and explored the evolution of 
either a core field threaded in the proton fluxoids, and a surface field buried by 
the diamagnetic accreting plasma. Here, we start with a population of NSs at 
the onset of the Roche Lobe Overflow (RLO) phase and include the physical 
evolution of a crustal field. The adopted parameters arc: 



Physical quantity 


Distribution 


Values 


Units 


NS period at t^^^ 


Flat 


1 


^ 100 


sec 


NS At at t^LO (a) 


Gaussian 


Log< Ho > 


= 28.50 ; (7=0.32 


G cm3 


m in RLO phase ^''^ 


Gaussian 


Log< m > 


= - 1.00 ; f7=0.50 


Me 


Minimum accreted mass 


One- value 




0.01 


Mo 


RLO accretion time ^''^ 


Flat in Log 


10^ 


, ^max (d) 

^ ^RLO 


year 


MSP phase time 


Flat in Log 


10^ 


3 X 10^ 


year 



(a) ^RLO ^ initial time of the RLO phase 

(''^ baryonic accretion rate during the RLO phase 

^'^^ a Maximum accreted Mass of O.SM© is permitted during the RLO phase 

max duration of the RLO phase; typical values: 5x10^ yr - 10^ yr - 5x10^ yr 



The values in this Table are derived partly from the observation of LMXBs and partly 
from the tracks in the /i — P plane calculated by Urpin, Geppert & Konenkov (1998) 
for a sample of pulsars, under different hypotheses for the binary evolution times, for 
the /z decay and for the angular momentum transfer at the magnetospheric radius. 

2. Results and Discussion 

To execute the synthesis calculation we built a Monte Carlo code, using typically 
3,000 particles. The statistical analysis is carried on only those NSs reaching 
the so-called "millisecond strip" at the end of recycling, the ones having period 
P < 10.0 ms and whichever value of the magnetic moment /x. In accordance 
with the values of Pmin and iimin (the weakest magnetic moment observed; 
Mmm = 7.3 X 10^^ G cm^), we divided our particles in four groups. Those filling 
the first quadrant in the millisecond strip (P > Pmin and f^ > jimin) behave as 
the known MSPs. Also the objects belonging to the second quadrant (P < Pmin 
and H > Hmin) should shines as pulsars (sec Burderi &; D'Amico 1997). The 
observability of the objects in the third quadrant (P < Pmin and /j, < jimin) as 
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Figure 1. Left: Statistical properties of the synthesized populations 
in the l-i — P plane, for < m >= 0.1 and rjj'^Q = 5 x 10^ yr. Full dots 
represent the sample of detected MSPs, while open squares represent 
the NSs in LMXBs (White & Zhang 1997). Right: Calculated distri- 
butions (with and without propeller) of millisecond NSs as a function 
of the spin period P, while covers the whole range of values. 

radio sources represents instead a challenge to the modern pulsar surveys. Most 
of these NSs will be above the Chen and Ruderman (1993) "death-line", and 
might have a bolometric luminosity comparable to that of the known MSPs. 
Thereafter we shortly refer as su6-MSPs to all the objects having P < Pmin 
and n above the "death-line". Objects in the fourth quadrant (P > Pmin and 
M < fJ-min) are probably radio quiet neutron stars (RQNSs), because they tend 
to be closer to the theoretical "death-line" , and they are in a period range which 
was already searched with good sensitivity by the radio surveys. 

The left panel of Figure 1 illustrates how the simulated NSs fill the ji — P 
plane, while the right panel shows how the EoS affects the shape of the period 
distribution. Figure 2 reports on the typical relative abundances (left) and gives 
indication on the role played by the propeller (right) in modifying the results. 

3. Conclusions 

The population synthesis calculation leads to a number of interesting predictions: 

A. It shows the presence of a tail in the period distribution of the simulated 
NS populations at periods shorter than 1.558 ms. 

B. The fraction of s?i&-MSPs estimated over the entire MSPs population 
varies between to ~ 50%. The EoS and the parameters for the recycling play 
both an important role in determining these percentages. 

C. For a mild-soft EoS and irrespective to the boundary condition at the 
crust-core interface for the magnetic field evolution, the recycling in LMBs gives 
rise to a NS distribution increasing toward short periods and a clear barrier is 
present at the minimum period for mass-shedding. 

D. For a stiff EoS, the distribution is flatter and displays a broad maximum. 
The maximum is located where most of the spin periods of the NSs in LMXBs 
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Figure 2. Distributions of the synthesized NSs, derived normaUzing 
the sample to the total number of stars with P < 10 ms. We have 
divided the H — P plane in four regions. The upper left number in 
each cross gives the percentage of objects having P < Pmin and ^ > 
IJ"min- Left: Distribution for the parameters of Fig. 1-left. Right: 
Distributions when a strong propeller is applied, as in Fig. 1-right. 

are found, as indirectly inferred form the quasi periodic oscillations seen in the 
timing analysis of their X-ray light curves. 

E. If NSs at the end of persistent accretion in a LMB experience a phase 
of smooth decline of the accretion rate, the magnetospheric propeller at the end 
of the RLO phase produces a depiction of fastly spinning NSs but, at least for 
the soft EoS, it preserves a distribution that peaks at periods ~ 1.5 ms. 

F. The models for the decay of a crustal magnetic field predict the existence 
of spun up NSs with very low magnetic moment: their period distribution is a 
neat signature for the physics at the crust-core interface. 

This investigation shows that the detection of sm6-MSPs represents a serious 
challenge to modern radio searches. 
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